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Abstract. We consider integrable discretizations of some soliton equations associated with 
the motions of plane curves: the Wadati-Konno-Ichikawa elastic beam equation, the complex 
Dym equation, and the short pulse equation. They are related to the modified KdV or the 
sine-Gordon equations by the hodograph transformations. Based on the observation that the 
hodograph transformations are regarded as the Euler-Lagrange transformations of the curve 
motions, we construct the discrete analogues of the hodograph transformations, which yield 
integrable discretizations of those soliton equations. 
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1. Introduction 

The study of discrete integrable systems has received considerable attention in the past decade 
(see, for example, HI). Ablowitz and Ladik proposed a method of integrable discretizations 
of soliton equations, including the nonlinear Schrodinger equation and the modified KdV 
(mKdV) equation, based on the Lax form [|2H31, and Hirota discretized various soliton 
equations such as the KdV, the mKdV, and the sine-Gordon equations by using the bilinear 
formalism (SHU. Following the pioneering work of Ablowitz-Ladik and Hirota, Date, Jimbo 
and Miwa developed a unified algebraic approach from the point of view of the KP theory 
[flQ] - [T6l . For other approaches to the discrete integrable systems, see, for example, H17H181I . 

It is known that there is a class of soliton equations which admits loop, cusp, and peak 
soliton solutions. Among them, some soliton equations, e.g. the Wadati-Konno-Ichikawa 
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(WKI) elastic beam equation, the Camassa-Holm equation, the Degasperis-Procesi equation, 
the short pulse equation and the Hunter-Saxton equation, appear as model equations of 
various physical phenomena lfT9H2~6ll . It should be noted that those equations are transformed 
to certain soliton equations which admit smooth soliton solutions through hodograph 
transformations li27l - [35ll . For example, the WKI elastic beam equation is transformed to the 
potential mKdV equation GTTl . 

Integrable discretization of those soliton equations has been regarded as a difficult 
problem until recently. Especially, a systematic treatment of hodograph transformations in 
discretizations has been unknown. Recently, some of the authors succeeded in integrable 
discretization of some equations in the above class of soliton equations by using the bilinear 
method, and it was confirmed that those integrable discrete equations work effectively on 
numerical computations of the above class of soliton equations as self-adaptive mesh schemes 
ll36l - [39ll . However, the method employed there is rather technical, so it is not easy to extract 
a fundamental structure of discretizations to apply this method to a broader class of nonlinear 
wave equations. 

On the other hand, the discrete integrable systems have been applied to discretizations 
of curves and surfaces, and this area has been recently studied actively under the name of the 
discrete differential geometry H401I411 . In particular, there have been intensive studies in topics 
related to curve geometry after the pioneering work of Lamb and Goldstein-Petrich l!42ll43ll . 
and then several frameworks for the motion of discrete curves have been proposed in various 
settings Il44l - l49ll . 

It is well known that the potential mKdV equation describes the motion of plane 
curves [|43ll . Recently, the authors considered continuous and discrete motion of discrete 
plane curves in the Euclidean plane and presented the explicit formula in terms of the r- 
f unction [|48ll49l . The hodograph transformation of the WKI elastic beam equation can be 
viewed as the Euler-Lagrange transformation of the motion of plane curves. From this fact, 
it may be possible to establish a discrete analogue of the hodograph transformation for the 
motion of discrete curves, and to discretize soliton equations having singularities naturally 
from a geometric point of view. 

In the present paper, we construct discrete analogues of hodograph transformations by 
the Euler-Lagrange transformations of the motion of discrete plane curves in the Euclidean 
plane. Based on them, we construct semi-discrete and fully discrete analogues of the WKI 
elastic beam equation, the complex Dym equation, and the short pulse equation. In Section 
2.1, we discuss the motion of plane curves described by the potential mKdV equation and 
the hodograph transformations for the WKI elastic beam equation and the complex Dym 
equation. We also discuss the relationship between the sine-Gordon equation and the short 
pulse equation. 

In Section 2.2, we introduce the discrete hodograph transformation for the continuous 
motion of discrete plane curves which are described by the semi-discrete potential mKdV 
equation. Then we construct the semi-discrete WKI elastic beam equation and the semi- 
discrete complex Dym equation. Using the same technique to the semi-discrete sine-Gordon 
equation, we construct the semi-discrete short pulse equation. 
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In Section 2.3, we consider the discrete motion of discrete plane curves which are 
described by the discrete potential mKdV equation. Then we construct the discrete WKI 
elastic beam equation and the discrete complex Dym equation by using the discrete hodograph 
transformations. In a similar way, we construct the discrete short pulse equation from the 
discrete sine-Gordon equation. 

In Section 3, we present a r-f unction which gives soli ton and breather solutions for these 
equations, and Hirota-type bilinear equations for the r-function. Section 4 is devoted to the 
conclusion. 



2. Motion of plane curves and hodograph transformations 

2.1. Motion of smooth curves 



Let y(s) be an arc-length parametrized curve in Euclidean plane R 2 . Then the tangent vector 
dy 



tt- satisfies 

os 



ds 



1. 



Thus y. admits the parametrization 



dy 
Ts 



cosO 
sin# 



(2.1) 



(2.2) 



dy 



The function 6 = 6(s) is called the angle function of y which denotes the angle of j- measured 
from the x-axis. We define the normal vector N by 



N 



-1 

1 

and introduce the Frenet frame 



dy 
Ts 



-sin# 
cos# 



F = (T,N), T=-^, 
os 

which is the orthonormal basis attached to the curve. The Frenet equation is given by 



d 


f ° 


—K 


—F = F 






ds 


K 






(2.3) 



(2.4) 



(2.5) 



where the function k = || is the curvature of y. The angle function 6 is also referred to as the 
potential function. Let us consider the following isoperimetric motion in time t: 



1f = f 

dt 







Kss "T 9 



"*« - T 



In terms of ^, (12.51 ) and (12.61 ) can be expressed as 



d 2 y 
ds 2 


-K 
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Ts' 
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K ss 


K 3 ' 


dy 
ds 



(2.6) 

(2.7) 
(2.8) 
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respectively. Then the compatibility condition of (12.51) and (12.61 ), or (12.71) and (12.81 ) yields the 
mKdV equation for k = k(s, t) H421I43H 

3 

K, + -K 2 K S + K sss = 0, (2.9) 

or the potential mKdV equation for 9 = 6(s, t): 

9, + 2^) 3 + 0*ss = 0. (2.10) 

The mKdV equation can be viewed as the governing equation of the Lagrangian description 
for the motion of the curves y in terms of the arc-length parameter s. Let us consider the 
Eulerian description of the same motion of the curves. To this end, we introduce the Eulerian 
coordinates 



y(s, i) 



x(s, t) 
v(s, t) 



f 

Jo 



cos 6(s',t) 
sin 6(s',t) 



ds' + 






(2.11) 



(2.12) 



and change the independent variables (s, t) to 

(x, t') = I cos 6(s', t) ds' + x , t 

For simplicity we write f as t without causing confusion. Let us write down the equation for 
v in terms of x and t. It can be easily shown that 



s(x, t) = 



N = 




1 + v\ dx, k(x, t) = 



(1+v^ 



V^I 



-V.v 
1 



T = 



1 



^R 



1 



Noticing 



it follows that 



- r = -,,N-^T, 



- k, = y t ■ N 



v t 



Vr 



(2.13) 

(2.14) 

(2.15) 
(2.16) 



+ vi 



by taking the inner product with ./V on both sides of (12.151) . By using j- = yTTvf, we see 
that 



v t 



Thus we derive 



-^V^i— <■■ 



v t = 



(1+Vl)2 



(2.17) 



(2.18) 



Introducing u = v x , we obtain the WKI (Wadati-Konno-Ichikawa) elastic beam equation lfT9T 

(2.19) 



u, = 



u x 



(l + u 2 )i 
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Therefore, (12.181) or (12.191) can be viewed as the governing equation of the Eulerian description 
for the curve motions given by (|2.5I) and (12.61) . 

We note that (12.111) is the hodograph transformation between the potential mKdV 
equation (12.101) and the WKI elastic beam equation (12.181) found by Ishimori [27J. The above 
discussion shows that the hodograph transformation arises naturally as the transformation 
between the Lagrangian and Eulerian descriptions from the point of view of geometry of 
plane curves. 

There is another equation related to the plane curve motions which is known as the 
complex Dym equation (this is often called the complex Harry Dym equation) H19U43II50I4551 . 
Introducing the complex variables r and z, we consider the transformation 

r(s,t) = e' fZIeis ' t \ (2.20) 

z(s, t)= f e^ 19 ^ ds' + zq, t' = t. (2.21) 



Then the potential mKdV equation (12.101) is transformed to the complex Dym equation H52U531 

r, = r 3 r zzz . (2.22) 

Here we set f = t without causing confusion. The geometric meaning of the complex Dym 
equation may be described as follows. The variables r and z are expressed as 

r = cos 6(s, t) + V^T sin 6(s, t) = x s + V^Tv,, z = x + V^Tv + z . (2.23) 

Identifying the Euclidean plane R 2 as C, we see that z is the position vector of the curve y, and 
r corresponds to the tangent vector 'j- s ■ By treating the Eulerian coordinates of y as a complex 
variable z, the complex Dym equation is nothing but the governing equation to describe the 
motion of the tangent vector |j. 

It is well known that the sine-Gordon equation 

6 ys = 4sin9, (2.24) 

belongs to the the same hierarchy as the mKdV equation H56U57I and that it describes a certain 
motion of plane curves [58J. It is possible to derive the governing equation of curve motion 
in the Eulerian description in a similar manner to the case of the mKdV equation. In fact, 
applying the transformations 



HI 



(*,/)= cos6(s',y)ds' + x ,y), (2.25) 



-f 

Jo 



sm9(s',y) ds' + v , (2.26) 



we obtain the short pulse equation H25U59J4621 

v.w = 4v + ^(v 3 )^, (2.27) 

where we set y' = y for simplicity. Again, we note that the short pulse equation (|2.27l) 
describes the same curve motions as the sine-Gordon equation by using the Eulerian 
description. The transformation (|2.25l) gives the hodograph transformation between them [33 , 
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2.2. Continuous motion of discrete curves 

In this subsection we discuss the semi-discrete equations arising from the continuous motion 
of discrete plane curves and the hodograph transformations among them. A map y : Z — » 
IR 2 ; 1 1-> 7/ is said to be a discrete curve of segment length ai if 

(2.28) 



(2.29) 





7/+i - 7/ 




= 1. 






ai 




We introduce the angle function if/ t of a discrete curve y by 


7/+i - 7/ 


COSl/f/ 




a, 


sin i/// 




A discrete curve y satisfies 


7/+i - 7/ P / n 7/ - 7/-i 


for Ki = if/i - if/i-i, where R(ki) denotes the rotation matrix given by 


^ /cos* -sinM 
\ sin a:/ cos/c/ / 


We set a/ = e(> 0), and consider the following motion of discrete curves 


dy, 1 
rf£ cos 


KI 

2 


*(-!) 


7/+i - 7/ 

e 



(2.30) 



(2.31) 



(2.32) 



Then from the isoperimetric condition (12.281) and the compatibility condition of (12.301) and 
(12.321) , it follows that there exists a potential function 8i characterized by 



<A/ 



@i+\ + 0/ 



K, = 



Qi+i - Q, 



7-1 



(2.33) 



2 ' ' 2 

and that 0/ satisfies the semi-discrete potential mKdV equation 

dBi 2 /0 /+ i-0 w \ 

-— = - tan . 

rf£ e V 4 / 

We note that K t = - tan | satisfies the semi-discrete mKdV equation 

f=!(i + £?)«..-*.>. 

It is possible to consider the Eulerian description of the curve motion defined by (12.301 ) 
and (|2.32l) . Noticing (12.291) and (|2.33l) . we introduce the Eulerian coordinates 

rin(^) 



(2.34) 



(2.35) 



7/GD = 






i=0 



6 COS I 



6 sin | 



+ 



X 
v 



(2.36) 



Then from (1Z291 . (12341) and (12361) . one can derive 

d_ s __ v/+i -vi I A/+i - A/ _ A/ - A/_i 
^ '" e \l+A /+1 A/ + 1+A ; A H 

d Si I A /+ i - A/ A/ - A/_i 

-TZ (V/+1 - V/) = — — + — 

dt, e\l + A, +1 A, 1+A/Am 



(2.37) 
(2.38) 
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where 

81 = X M - X u A, = sgn (v, + i - Vi) 



V/+1 - Vi 



e + 6, 
since e + 6/ = e(l + cos ifri) > 0. Note that v/ and 61 satisfy 



V/+1 - vi 
e + Si 



V/+1 - Vi 



5, 



+ 1-1 - 1. 



From (12371) and (12381) . we obtain 
d I v /+ i - vi 



dt\X M -X,) e{ \X M -X, 



2 \ 



A /+1 - A, A, - A M 



(2.39) 



(2.40) 



(2.41) 



1+A /+1 A, 1 + A/Am, 
The system of (|2.37l) . (12.381) and (12.391) is nothing but the semi-discrete WKI elastic beam 
equation. We remark that (12.361) can be regarded as the hodograph transformation between 
(12.341) and the semi-discrete WKI elastic beam equation. Note that the angle function 



^£& satisfies 



cosi/// = 



Xi+i - Xi 



sini/// = 



v/+i - v, 



tani/// = 



V/+1 ~ Vj 

Xi+i - X{ 



Thus (12.411) can be rewritten as 

_fc = -(tan-^— +tan— ^— 

Equation (|2.43l) with the discrete hodograph transformation 



(2.42) 



(2.43) 



i-\ 



i- 1 



X,(0 = ^ e cos ifrj(0 + X , v,(0 = Yj e sin <M£> + v " ' ( 2 - 44 ) 

7=0 7=0 

can be also regarded as the semi-discrete WKI elastic beam equation. In the continuous limit 



e -> with s = el + £ and t = ~£, (12.431) and (12.441) converge to 



e t + -(e s ) 3 + e sss = o , 



and 



x(s, t) 



-f 

Jo 



cos 8(s', t)ds' + x , v(s 



,t)= f 
Jo 



sin 6(s', t)ds' + v ■ 



(2.45) 



(2.46) 



which give the (potential) WKI elastic beam equation (12.181) (see Appendix). 

One can construct a semi-discrete version of the complex Dym equation as follows. In 
view of (12.231) . (12.331) and (12.361) . it is natural to introduce the complex variables r/(£) and 
Z/(£)by 



/-i 



n(0 = e^' = e ^ ^ , ZU) = X£0 + V=lvi(0 = £ « ^ 

7=0 



+ Z n 



(2.47) 



(2.48) 



Then we have from (12341) and (12471) 

dn _n ( n+\ - n n - r^ \ z /+ i - Z/ _ 

d£ e\r/ + i+r z n + ri-J' e 

which is the semi-discrete complex Dym equation. The geometric meaning of (12.481) can be 
described as follows: under the identification of M 2 as C, Z/ is the position vector of the curve 
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8 



ji, and rj corresponds to the segment vector y/+ '~ y ' . Then (12.481 ) is the governing equation 
describing the motion of the segment vector in the Eulerian coordinates of y/. To take the 
continuous limit, we use the angle function <///. Then the semi-discrete complex Dym equation 
(12.481) is rewritten as 



d 1 / iff M -ij/i </// - if/i-i 

— iff, = - tan + tan 

dC e\ 2 2 



(2.49) 



and 



i- 1 



ZU) = £ ee ***#> + Z , rU) = e +**<> 

7=0 



(2.50) 



In the continuous limit £ — > with s = el + t, and t = -t-^, (12.491 ) and (12.501) converge to 



e t + -(o s f + e sss = o, 



and 



Z(s,t)= [ e ^^'"W +zo, r(s,t) = e ^ e{s ' t: 
Jo 



(2.51) 



(2.52) 



which give the complex Dym equation (12.221) (see Appendix). 

Now we construct the semi-discrete short pulse equation. To this end, we consider the 
semi-discrete sine-Gordon equation 



-t(0m ~ Si) = 46 sin I 



(2.53) 



Similar to the continuous case, the semi-discrete sine-Gordon equation (12.531) can be regarded 
as describing a certain motion of discrete plane curves. Therefore, we may expect that the 
application of the same transformation as the case of the semi-discrete WKI equation to the 
semi-discrete sine-Gordon equation (12.531) yields the semi-discrete analogue of the short pulse 
equation. In fact, by using the transformation 



yiiy) = 



Xi(y) 

viiy) 



■Z 

7=0 



6 COS | 



esinl ' 2 ' 1 



+ 



X 

V() 



we obtain the semi-discrete short pulse equation 

j{X M -X 1 ) = -2{v 2 M -v 2 l ), 

d 

-tOm-i - vi) = 2(X M - X/)(v,+i + v/). 
dy 

We note that the following relation also holds from (12.541) 

* 2 <X M -X^ 2 



/ v/ + i -v/ r / x M -xrf 



= l. 



6 / \ £ 

From (|2.55l) and (12.561) . we obtain 



d I v /+ i - vi 



dy\X M -X, 



2(v /+ i + v/) + 2 1 -^ — -^ | (vi+i + vi) . 



(2.54) 

(2.55) 
(2.56) 

(2.57) 

(2.58) 
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In order to take the continuous limit, we assume the boundary condition X; = v/ = for / < 0, 
which is consistent with (|2.54l) . Then the continuous limit e — > (i.e., X t+l - X t — » 0) gives 
V/+i - vi dv v /+ i + V[ 

*dx> ~^^ v > 



Xi+i - Xi 

dX i = dXo 
dy dy 



+ y d(X j+l - Xj) _ dX { 
7=0 







dy 



dy 



./=<> 



d d dx d d 2 d d n 2 ® 

dy dy' dy dx dy' l dx dy' dx 



Consequently, (12.581 ) converges to 



{d y - 2v 2 d x )v x = 4v + Avv 2 x > 



which is nothing but the short pulse equation (I2.27I ). 



2.3. Discrete motion of discrete curves 



9X r, 2 

— = -2v 
dy 



(2.59) 



Now let us recall the following discrete motion of discrete plane curve y% introduced by 
Matsuura [|49l : 



r, 



;/+l 



Yn 



a n 


y m - 
Yn+\ 


y: 


a„ 




jT 1 - 


Tn 



= 1, 

y m _ m 

•■ RiO - — ^±, 

®n-\ 

m _ m 
I n+\ in 



R«) 



(2.60) 
(2.61) 
(2.62) 



"m (*n 

where a„ and b m are arbitrary functions in n and m, respectively. Compatibility of the system 
(12.601) — (12.621) implies the existence of the potential function ff" defined by 



ffn+l an an _ an 

n _ n u n +[ m _ U n+l U n-\ 

UJ m — _ , K„ — 



2 " 2 

and it follows that 9% satisfies the discrete potential mKdV equation [[68 



( ff^-ff^ b + a 19^-9? 
tan | = tan | - 



n+\ 



4 / b m -a n \ 4 
Note that the functions t//™ and (p™ can be expressed as 



€ 



an , an 
6 n+\ + 6 n 



<f>n 



an+i _i_ an 
m _ U n + V n 



(2.63) 



(2.64) 



(2.65) 



Note also 



and 



yn _ yn 
i n+1 in 

a n 



cos i//'" 
_ sin if,™ 



-.m+l _ yn 
in In 



r n = 



in 
V 



n-\ 

=z 

7=0 



vm vm 

A 7+l " X j 

m m 

V , — V . 

7+1 7 



+ 



X" 



COS 0™ 

sin0™ 



rc-1 



=Tj a j 



7=0 



COS I//'" 

sini//" 1 



+ 



X" 



(2.66) 



.(2.67) 
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From the discrete potential mKdV equation (12.641) and the hodograph transformation (12.671) , 



we obtain 



where 



A' 



m+l 



h + a Y m - A m 

u m -r u n ± n is- n 



l + A™ +1 r;f b m - a n l + r™A;? 



K = sgn(v;r +1 - o 



r^sgnCv^-v^ 1 ) 



*n+l K n 



On + (*£, ~ X-) 

v m+l _ y m 



"n+l y « 



..m+l .jn 



(2.68) 

(2.69) 
(2.70) 



Z? m + (Z» +I - X™) 

since a„ + (X^ +1 - X") = a„(l + cos C) > and b m + (X;" +1 - X") = & w (l + cos O > 0. We 
note that v™ and X'" satisfy the following relations 



-ym vn 

A n+l ~ A n 



a n 



X 



m+l 



X" 



+ 



+ 



n+l 



a n 



= 1, 



= 1. 



(2.71) 



(2.72) 



b \ h 

u m i \ u m 

To construct an explicit form of the discrete WKI elastic beam equation, we consider an 
identity 



Substituting 



gV^r^v^Tc _ e v=T^ lg v=r« 



e V^Tc =cos ^+ V^Tsin^ 



(2.73) 



e^* = cosC + V-Fsin0™ = -*■ 



ym ym ,,m ..m 

A n+1 ~ A n ^T^- V »+l ~ V n 

ym+1 ym 



+ 



a n 

-..m+l ,.m 



b b 

into (12.731) . we obtain the system of two discrete equations from the real and imaginary parts, 
respectively, which should be considered together with the constraints (12.711 ) and (12.721) . To 
incorporate (12.711) and (|2.72l) . an easy way is to consider the identity 

e V=T#r V V=T* = e V=TCi g - V=T« f (2 .74) 

instead of (12.731) . We then obtain the following system of two discrete equations 



/ym+1 -ym+l\f-ym+\ -ym\ . /-m+l ,,m+l\/,,m+l ,,m\ 
^ A n+1 _A n )\ A n ~ JL n) + \ V n+l ~ V n A V n ~ V n> 



rm+\ 



_ /ym+l 
~ ( A n+1 



K + i)(K 



+i 



K) + (d 1 



CiXCi-O. 



V A «+1 _A n AVn _V n^ _( , A n " A Jl V »+l _ V n -> 
/■m+l -i^m \/,.m ,.m\ /y^n ymw. .m+l 



= (ez - x^iXCi - c) - «i - ^rxv 



n+l 



"n+l 



), 



(2.75) 
(2.76) 



which is simplified to 



ra+l 

m+l 



V A «+1 " A 

+ (v n+ i 

/ym+1 . ym ym+1 

V A „+1 + A n+1 _ A n 



m+l\ 



Xm+l . yfflvyffl yffl+h 

„ + A „ A A n+l ~ A n ) 



Ci - v 



m+l 



n+l 

m 



+v:xci-c i ) = o, 



ym\/ m+l 
A n )\ V n+\ 



v;+i + v 



,m+ 1 



v m ) 



/ym+l ym . ym+1 \7fn\/.m+l , .m ,/w+l ,,m\ _ n 

_ 1 A »+1 -A n+1 +A n ~ A nA V n+l + V n+\ ~ V n ~ V n ) ~ U 



(2.77) 
(2.78) 
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Note that the second equation (12.781 ) is further simplified to 

(X£i - XT'XCi - v «) - (C - *?XCi - vT 1 ) = ■ (2.79) 

We remark that the system (12.771) and (12.781) can be solved explicitly in terms of V^\ and 
X m +} in the form of rational functions of v"\ v n \ „ v"?, 1 , X"', X"' and X' n+1 . Therefore the 
system (12.771) and (12.781) (or (12.791 )) can be regarded as the explicit form of the discrete WKI 
elastic beam equation. Note that we can obtain (|2.68l) . (12.691) and (12.701) by replacing X™ by 
X™ = X™ + 2"ld aj + Z'jS/ bj in (12/751) and (12J61) and then dividing (12J61) by (1X751) . By using 
the potential function 0™, the discrete WKI elastic beam equation can be written as 

, an /^U^ t R"'-C.\ (2 . 80) 



\ " ' / "m ®n \ 



n- 1 f am , an 

7+1 ' ' '+X™. (2.81) 



X- = ^a ; cos - 

7=0 ^ 



Setting 

a + b a — b 

g = (n + m)£ , I = n-m, a n = a , b m = b , 5 = — - — , 6 = — - — , (2.82) 

and taking the continuous limit 6 — » 0, (12.801) and (12.811) are reduced to 

^ = ; tan (— M' (2 - 83) 

X,(£) = J, a J cos 2 + W) > (2-84) 

which are transformed to the semi-discrete WKI elastic beam equation (12.371 ), (12.381) and 
(12391) . 

Let us consider a discrete analogue of the complex Dym equation. Introducing 



»"', , +C 



c = gV =T-V L , (2.85) 

n-\ n-\ 

Z: = X>: + V=Tv? = ^flycos^ + V=T2>sin^ +Z- 

7=0 7=0 

n-1 

= 2«^+Zo. (2-86) 

7=0 

and using the discrete potential mKdV equation (12.641) . we derive the discrete analogue of the 
complex Dym equation 



m+l „m u i „ ,jn vJn »jn+l rfn 

(2.87) 



r n Pn °m + "n Pn r n ' n l^n+l 



t 4n+\ _i_ f^m U n f^m _i_ v m ' .yn r m ' 

'n T Fn ^m "« rn "*" 'n rn ' n 

Z^-Z: = a n C (2.88) 

where p™ is an auxiliary variable defined by 

p m n = e^^r^ . (2.89) 
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Similar to the semi-discrete case, the system of (12.871) and (12.881 ) describes the motion of 
segment vector of the curve y™ in the Eulerian coordinates of y™. It should be noted that by 
introducing Q'" by 



q: 



.V^Te;; 1 



we have an alternate form of the discrete complex Dym equation 

r m = lo m O m 



QZ1+J& bm ~ an ^2^+ VS"' 

Z' n = aJ! . 



Using ff£, the discrete complex Dym equation can be written as 



tan 



am+l an 



b m + a n (ffn + ~0 n n \i 
— — — tan 

"m "n \ T 



n-\ 



7=0 






+ Z" 



Setting 



a + b 



£ = in + m)6 , I = n - m, a n = a , b m = b , 8 = 

z 

and taking the continuous limit c> — > 0, (12.941 ) and (12.951 ) become 

dQi 2 (Bi+i - 61-1 
-— = - tan 

d£ e \ 4 



i-\ 



Zi(0 = J] aje 



An - «; + l(f)+«/f) 



+z (O, 



(2.90) 

(2.91) 

(2.92) 
(2.93) 

(2.94) 
(2.95) 

(2.96) 

(2.97) 
(2.98) 



./=(> 



which are transformed to the semi-discrete complex Dym equation (12.481 ). 

We next construct the discrete short pulse equation. Consider the following discrete 
motion of plane discrete curve y k n [48J 



1 n+1 'n 



a,, 



= 1, 



k k k k 

^tj+1 — in nt ,k\ in ~ 'n-\ 



a„ 



a n -\ 



In ~ * Tn = —^ R (-<K) ■ 

Cfc a n 



1 
-1 



(2.99) 
(2.100) 
(2.101) 



where a n and c* are arbitrary functions in n and k, respectively. Compatibility of this system 
implies the existence of the potential function 6^ defined by 



rk _ rk 
k _ n+1 U n-\ 



k _ ^ +1 + C 



(2.102) 
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and it follows that 8 k satisfies the discrete sine-Gordon equation 

. / C-C-^ + ^ \ ■ (C\ + Ci + ^ +1 + # 

sin = a n Ck sin 

Note that the functions ifj k n and (p k can be expressed as 

/k c + g * e 1 - # 

¥n = o ' ^n = o ■ 
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(2.103) 



(2.104) 



Note also 



and 






cosi// k 
sin^ 



rr-Sy k n 






COS^ 

sin<^ 



rf, 



V* 



0-1 

7=0 






+ 


[41 


«-i 


COS (/>y 

sini//^ 


+ 


r y* i 

v k 



(2.105) 



.(2.106) 



From the discrete sine-Gordon equation (12.1031) and the hodograph transformation (12.1061) . 
we obtain 



A* +1 - A* = a n c k (Y k ; 1 + r*) , 



(2.107) 



where 



A* = sgn(v* +1 - v*) 



V* , — V* 

"n+1 K n 



r^sgn^+v^) 



a» + (X* +1 -X*) 



*+l + i 



v k , -v k 



i + W - x k ) 



a 



a n + (X k +l - X k ) 

V„ +1 +Vn 

i + (*r - 3D ' 



(2.108) 



(2.109) 



since a„ + (X* +1 - x£) = a„(l + cos if/ k n ) > and ^ + (X k+l - X k ) = ±(1 + cos ^) > 0. We note 
that v k n and X* satisfy the following relations 



>7, + .-M + /^ + i-« 



U/i 



< +1 +v k „ ) 

c k J 



a„ 



+ 



± 



k\ 



= 1. 



(2.110) 



(2.111) 



We now construct an explicit form of the discrete short pulse equation. Similar to the case of 
discrete WKI elastic beam equation, we consider the identity 

g V=T*£VV=Trf = e V=T^V^* + , , (2.112) 



with 



e^ 1 ^ = cos(//f, + V-Tsin^ = 



X n+ l ~ X n /— j- V n+l ~ V n 

a n 

»-rl + v k 



Cl n 

x k+l - xt 



>*** = cos^ + V^Tsin^ = ^ n ; ^ M + V^T ^ 






_l_ 
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We then obtain the following system of two discrete equations 



rk+l 



■k+l\ 



\ A n+l ~ A n+l ~ A n + A n)\ A n+l ~ A n > 



(x k : + \ 






- (v£l + V «+l - V 



jfc+1 



-v*xci + vr 1 ) = o, 



X n+l + X[ 



-k+\ 



*m:\ 



'n+1 



v k+l + v k ) 

y n y n> 



- (x k+ + \ + x k n+1 - x k+1 - x k M:\ + v* + i + < +l + 



0. 



Note that the second equation (12.1141) is further simplified to 

(X k +1 - X k :'){v k n \\ + v k n ) + (x£J - X>* +1 + v* +1 ) = • 
Replacing X k by X k + £*:' -±, (1XTT31 and (12TTT41 become 



TO-*; 



l n+l 



x! +1 - - 



- (K+i + K + i ~ v 



c k 

i+i 



— +X: 



n+1 xr+X*)\X; 

k+l ■*■* ~ k+l -v k n M +l + < +l ) = o, 

k 



k+l 

71+1 



V^ _i_ vk+\ 
A n+l + A n 



x k M:\ 



V 



11+ 1 



y n y n> 



- (xjri + 4.1 - *t +I - *£)0fti + v£.i + v^ +1 + v k n ) 



o. 
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(2.113) 
(2.114) 

(2.115) 



(2.116) 



(2.117) 



Note that this form was obtained in [37| by using the bilinear method. Taking the continuous 
limit Ck — » of (|2.116l) and (|2.117l) . we obtain the semi-discrete short pulse equation (|2.55l) 
and ([236]). 



3. r-function and soliton type solutions 

In this section, we list the r-function and the bilinear equations which give rise to the 
soliton and breather type solutions to the equations and curve motions discussed in Section 2. 
Although they have been already discussed in |[37l|48l[69]|, we c °U ect an d present the results 
for completeness and the convenience of readers. It should be remarked that all the solutions 
can be expressed in terms of one r-function. 

The solutions can be expressed in the following form: 



0%(s,t,(,y) = 



Y H \(s,t,£,y) = 



m.k 



V^T 



log- 



m,k 
~n,l 



t m,k 
n.l 



\ n,l I v 



Here, the r-function r n ) (s, t, £, y) is given by [48 J: 



m,k/ 



r n )(s,t,t;,y) = exp 



n—\ m—\ k—\ -. \ " 



e'" + e^ 



(3.1) 



(3.2) 



(3.3) 
(3.4) 
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n—\ m-\ k—\ i \-\ 



n—i m—i k— i / v — i 3 Pi I 

e* = a iP \il - e Pi y> f](l - a^pd- 1 ]> ~ h ^ W V " ~ e"** 1 "^'**', 

n' m! k' * "' ' 

«-l m-1 fc-1 / v -1 3 Pi 

^ = A(-p/) ; (i + epo'' ri (i + an,piyi ri (i + bm ' pi) ~ i n ! + — e ~ /,v+ ^"^ 



The parameters are chosen as follows: 
(i) N-soliton solution: 

/?,-, a,- 6 R, #■ 6 V^Tm (i =1,...,N). (3.6) 



(3.7) 



(ii) M-breather solution: 

N = 2M, p h a u ^ 6 C (i = 1, . . . , 2M), 

P2j = P* 2j -i, a 2j = a 2j _ i , fo j = -p 2J _ i ( j = 1 , . . . , M) . 

The bilinear equations which are necessary to recover the equations and curve motions are 
given as follows. Note that we only show the relevant independent variables, and other 
variables can be regarded as parameters. 
Continuous case [f37l|48]|: r = t(s, t,y), 

^D s D y r-T=-(T*) 2 , (3.8) 

D] t ■ t = 0, (3.9) 

(D*+D t ) t-t* = 0. (3.10) 

Semi-discrete case |J371|69): T = T /(£»)0» 

D^Xi-t] = — (t^tj+i -t; +1 t/_i), (3.11) 

2 D ( D y t '-ti = -r*M T U> ( 3 - 13 ) 

(3.14) 





D y T/ + i -T/ = -eT /+1 T,. 


iscr 


ete case [[3711481: r = rf'^v). 




D y T n+l ■ r n = -a n T n+l r n , 




p> m+l,k m,k _ 1 *m,k *m,k 




r) m,k+l *m,k _ *m,k+\^m,k 
u y T n ' ' n ~ „ n T n ' 



(3.15) 
(3.16) 

(3.17) 



Ck 



4. Conclusions 

In this paper, we have discretized several soliton equations which admit loop type soliton 
solutions through the discrete analogues of the hodograph transformations based on the 
geometry of plane curves. More concretely, we have constructed semi-discrete and fully 



(3.5) 
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discrete versions of the WKI elastic beam equation, the complex Dym equation and the short 
pulse equation, and presented the r-function which gives rise to the soliton and breather 
solutions. 

Geometric consideration is effective for discretization of soliton equations which admit 
soliton solutions with singularities, and it may be also applicable to other soliton equations. 
For example, the soliton equations arising from the curve motions in the Minkowski plane are 
of so-called "defocusing type" which have nonlinear terms with different signs compared to 
the equations discussed in this paper, and it is known that the structure and behaviour of the 
solutions are quite different. The geometric consideration may also be useful to discretize this 
class of equations. This problem will be reported in a forthcoming paper. 

Appendix: Hodograph transformations 

The WKI elastic beam equation: 

A conservation law of the potential mKdV equation (12.101) is given by 

(cos#) r + l-(#,.) 2 cos#-#, s sin#) =0. (A.l) 

Consider the hodograph transformation |[2~7Tl 

x(s, t) = I cos 0(V, t)ds' + x , t'(s, t) = t, (A.2) 

Jo 

which leads to 

d d d d I 1 , \ d 

— = cos#— , — = — + \e, s sme--(9 s ) 2 cose — . (A.3) 

ds dx dt df \ 2 V J ) dx 

Applying (IA.3I ) to (IA.1I ), we obtain 

^+cos 2 0(sin 9) xxx = 0, (A.4) 

which can be rewritten as 

(tan0) r +(sin0);ccc = O. (A.5) 

Introducing a new dependent variable v(s,t) = JC sm8(s',t)ds' + v (note tan# = v x and 
sin 9 = v x l yjl + (v t ) 2 ), (|A.5I) is transformed to 



/ 



Vfx + 



0, (A.6) 



V 1 + (^) 2 L, 

which is the (potential) WKI elastic beam equation (12.181) . 

The complex Dym equation: 

A conservation law of the potential mKdV equation (12.101) is given by 

{e+*\ + foe.fe**' + ^\9 ss e^A = . (A.7) 

Consider the hodograph transformation H52U53H 

z(s, t)= f e ^^'^ds' + zo , t'(s, t) = t, (A.8) 
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which leads to 

5-^5- JrH> 2 ^"-^ w '^%- (A ' 9) 

Applying (IA.9I ) to (IA.7I) and introducing a new dependent variable r = e^ e , we obtain the 
complex Dym equation 

r f +r 3 (r) zzz = 0. (A. 10) 

Note that 

z = x+ ^/-[v, r = — . (A.11) 

as 

The short pulse equation: 

A conservation law of the sine-Gordon equation (|2.24l) is given by 

(cos 0). + (-£-) =0. (A. 12) 



8 
Consider the hodograph transformation 03311631 - 16511 

x(s,y)= I cos 6(s', y)d s' + x , y'(s,y)=y, (A. 13) 

Jo 

which leads to 

A = C0S #A A = A_^!A (A 14) 

ds dx ' dy dy' 8 dx 

Introduce a new dependent variable 

r r s QyAs',y) i 

v(s,y)= sm6(s ,y)ds' + v = -—— ds + v = -G y , (A. 15) 

Jo Jo 4 4 ' 

then it follows 

v x = tan . 
Applying (|A.14t to (1A.12K we obtain 

2v 2 — )cos6> = -4vv x cos6. (A. 16) 

o/ ox J 

this can be rewritten as 

#7-2v 2 |-)^- = 8vv,^-. (A.17) 

ay' ox) cos 1 6 cos 2 



From 

1 



, = 1 + tan 2 9 = 1 + v 2 
cos 2 



it follows that 



W ~ 2v2 ~oh) (l + v * } = 8vVv(1 + v%x) ' (A - 18) 



which is nothing but the short pulse equation 

3 v /xx 



4v + -](v 3 ) . (A.19) 

A \ I xx 
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